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Abstract. Quantum feedback control protocols can improve the operation of 
quantum devices. Here we examine the performance of a purification protocol when 
there are imperfections in the controls. The ideal feedback protocol produces an x 
eigenstate from a mixed state in the minimum time, and is known as rapid state 
preparation. The imperfections we examine include time delays in the feedback loop, 
finite strength feedback, calibration errors, and inefficient detection. We analyse 
these imperfections using the Wiseman-Milburn feedback master equation and related 
formalism. We find that the protocol is most sensitive to time delays in the feedback 
loop. For systems with slow dynamics, however, our analysis suggests that inefficient 
detection would be the bigger problem. We also show how system imperfections, such 
as dephasing and damping, can be included in model via the feedback master equation. 
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1. Motivation 

It is well known that imperfections in feedback control protocols — such as finite strength 
feedback, inefficient detection, and loop delays — severely degrade the controllability of 
some systems [H EJ |3j E]. In this paper we are interested in these types of control 
imperfections, as well as imperfect operations and non-negligible damping. In the 
interest of keeping this analysis manageable we restrict ourselves to a single, very simple 
measurement model, which we now describe. 

Measurement-based quantum feedback protocols often consider a system (we will 
examine a qbit) that has an interaction with a probe (a current or field) that results 
in a continuous quantum non-demolition (QND) jl] measurement||| of the system. The 
change to the system state p during continuous measurement of an observable X is 

dp = dtC c [X] p = 2-fdtV [X] p + ^y2yqdW(t)n [X] p. (1) 

The first term, T> [A] p = ApA^ — ^(A^ Ap+pA^ A), describes the measurement backaction 
on the state. The final term, % [A] p = Ap + pA< — Tr [(A^ + A)p\ p, represents the 
refinement to the observer's state of knowledge due to the measurement process [SlH]- 
The positive coefficient r\ is the measurement efficiency An efficient measurement 
corresponds to 77 = 1, while an inefficient measurement corresponds to 77 < 1. The 
measurement result in the interval [t,t + dt) is dR{t) = yffir/dt(X) + dW{t) / \/fj, where 
dW(t) is the Wiener process appearing in Eq. ([T]) [7]. The integrated current is given 
by R(t) = f*dR(s). 

In this situation the most general feedback strategy consists of applying a 
conditional Hamiltonian (which may be a functional of the entire record up to that time) 
to the system to effect the desired outcome. An important class of feedback strategies 
are those where the functional of the record is also a function of the best estimate of 
the system state. These strategies are called Bayesian [8] or state-based feedback [9]. 
In general, implementing feedback involves engineering driven Hamiltonian terms of the 
form —idtak{t)[Hk, p] to Eq. (Op). It is common to assume that the controls available 
are sufficiently strong and fast that one may consider directly controlling the state or 
measurement basis through a unitary U. (Note that this does not break any rotating 
wave approximations as we only require \a\ 3> |7|.) In this case the conditional state 
after the measurement and feedback is pit + dt) = p + dtC c \X] p in the Heisenberg 
picture with respect to the control unitary, where X = W(t)XU(t). 

Another important class of feedback control strategies can be described by the 
Wiseman-Milburn Markovian feedback master equation (FBME) [TU 1 [TT | H]: 

> = £ H--S 

I A simplified definition of a QND measurement is one where the measured observable commutes with 
the time evolution operator of the combined probe field and system, i.e. [U, X] = 0. See Refs. [HE] for 
a though treatment. Here we do not consider the control unitary to be included in this U — its action 
(after the measurement) certainly does not commute with X, as we will see. 



V27 



(XF + FX),p 



+ 2 7 £> 



X 



iF 



Quantum feedback for rapid state preparation in the presence of control imperfections 3 

+ l -^-V [F] p. (2) 
V 

This can be derived by considering instantaneous feedback of the current via the 
feedback Hamiltonian H^(t) = FdR(t)/dt, and averaging over the measurement noise 
to obtain a deterministic equation [HJH]. Here F is a Hermitian operator, heuristically 
chosen to effect the desired outcome. Markovian feedback is not as powerful as Bayesian 
feedback; it does however have the advantage that it is experimentally and theoretically 
simpler. 

In this article we systematically examine the effect of imperfections in quantum 
feedback protocols with a specific example - rapid state preparation (which is very 
closely related to Jacobs' rapid purification of a qbit [121 H31 EH])- We have chosen this 
protocol for a number of reasons. Firstly, rapid state preparation (RSP) has all of the 
ingredients of more complicated control protocols and experience shows that analysis 
of RSP often admits analytical solutions. Secondly RSP, as we shall soon see, is 
a highly quantum protocol as it requires continually creating (or stabilizing) a state 
that is a superposition of the pointer basis [15] - what would, for a larger system, 
be called a Schrodinger cat state. Consequently studying the effects of imperfections 
on quantum control protocols will help delineating in what regimes control will be of 
practical importance. Finally, much of the excitement about quantum information 
theory stems from the speed-ups over classical protocols. The RSP protocol is one 
quantum control protocol that offers such a speed-up. 

In the ideal limit, the RSP protocol can be achieved with Markovian feedback. To 
maintain analytic results, in this article we work within the framework of Markovian 
feedback; if we had worked in the Bayesian feedback framework we would have been 
forced to rely on numerics to obtain any results. There is also the possibility of using 
risk sensitive or robust quantum control [16] to examine the problems, but such analysis 
is also likely to be more complicated than the analysis we present. 

This article is organized as follows. In Sec. |2] we present a protocol for state 
preparation that minimises the use of feedback control to the application of one 
conditional unitary. In Sec. |3] we present the general theory for RSP allowing for 
continuous modulation of a feedback Hamiltonian under ideal conditions. Then each of 
these imperfections described above are examined in turn in Sec. HI Finally we conclude 
by comparing the relative detrimental effects of each imperfection. 

2. State preparation with one conditional unitary 

We begin by presenting the simplest state preparation protocol for producing a qbit in an 
x eigenstate, based on continuous monitoring of the z component of angular momentum 
i.e. J z . (We choose to use J z rather than the Pauli operators as it will make comparing 
our results to other results in the literature easy. This is because much of the work in 
the literature is for qudits.) It is obvious that this measurement in isolation will only 
produce eigenstates of J z . Thus, to produce an x eigenstate at time tf we assume that 



Quantum feedback for rapid state preparation in the presence of control imperfections 4 



it is possible to perform an instantaneous conditional unitary at that time to rotate the 
Bloch vector to the desired x eigenstate. This is the obvious method for minimising the 
use of feedback control. We assume here (and in all subsequent protocols) that initially 
the qbit is in a completely mixed state. 

We quantify the effectiveness of this protocol by the probability that it would be 
verified (by a hypothetical projective measurement) to have produced an x-eigenstate 
at time tf, as motivated in Ref. [H]. This probability is linear in the x component of 
the Bloch vector after the final conditional unitary rotation. This in turn is equal to 
the length of the Bloch vector prior to that rotation. We take this average length as the 
figure-of-merit for all of our protocols. 

In the present case, it is possible to obtain a closed form expression for the average 
Bloch vector length as a function of the measurement record by using linear trajectory 
theory [TTJ [HI [19]. This involves using a linear, and necessarily unnormalized, version 
of Eq. (P, called the linear SME 

dp = dtC c [J z ] p = 2-fdtV [J z ] p + dR H [J g ] p, (3) 

where H[A] p = Ap + pA^ . Here the bar denotes the lack of normalization. Given that 
p(0) and our observable J z commute, the solution to the SME is [6] 

p(R, t) = exp(-4 7 J 2 z t) exp(2 y/¥fJ z R{t))l/2, (4) 

where R(t) is the integrated photocurrent. Omitting the time dependence on R for 
compactness, the solution for a qbit can be expressed as 



e 



-7* / e V2yR o 



= — ( o e -^R )• ( 5 ) 



The norm of Eq. (jl]) is J\f — Tr [p(R, t)] = e 7 ' cosh {^/2PyR). The normalised state 
is then p(R,t) = p(R,t)/J\f. The actual probability density for the result R is 
V(R, t) = Tr [p(R, t)} P(R, t). Here 

p -R 2 /(2t) 

P(R,t) = ^= r - (6) 

has been called the ostensible probability for the result R [H] . The ostensible probability 
is the probability with which the record R should be generated if Eq. ([3]) were to be 
used in place of the normalized SME (CQ). 

With these results we may now calculate the average Bloch vector length. First 
we note that the Bloch vector length is given by 2A max [p] — 1, where A max [p] is the 
larger eigenvalue of the normalized state p. From Eq. (JHJ), this will correspond to the 
probability of being in the +z eigenstate or the —z eigenstate, depending on whether 
R is positive or negative. Using the symmetry of V(R) and Eq. (jSJ) the maximum 
eigenvalue is 



POO 

(A max (t)> = 2 / A max [p(t)] Af(R)P(R)dR 
Jo 

e - {R -^tf/2t dR 



1 
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= i[l + erf(V7*)]. (7) 
When t 3> 7 _1 (the long-time limit) this expression can be approximated by 

(A m „(i))LT = 1 - ~, (8) 

where the subscript LT denotes that this expression is only valid in the long-time limit. 
We can relate this expression to (x m fb), the mean length of the Bloch vector (or its 
x-component after the final rotation) by 

(ZmfbW) = 2(A max (t)> - 1. (9) 

Here the subscript 'mfb' signifies that this protocol minimises the use of feedback control. 
We also refer to this as the open-loop strategy because the conditional unitary happens 
after all measurements so it is feedforward. 



3. Ideal rapid state preparation 

Consider again a qbit in which J z is continuously monitored but now we assume that 
for arbitrarily strong feedback-induced rotation around the y axis is available. This is 
all that is required for control given that the measured observable picks out the z axis, 
and that the aim is to maximize (x(tj)). As we will see, the optimal feedback protocol 
in this situation is of the form described by Eq. (J2J), where H^{t)dt = FdR(t), provided 
F(t) itself is allowed to be time dependent: F(t) = Q(t)J y . Equation (E]) will thus 
provide the basis for the following study. 

Given the above senario and allowing for a time dependent feedback strength, Q(t), 
and assuming efficient measurements (77 = 1) , equation |2] becomes 

' v/27 J 

From Eq. ( TTOj) we may find the the equations of motion for the Bloch components. This 
is achieved by taking the trace of Eq. ( flOl) with the component of interest, for example 
x = Tr [a x p\. The resulting equations are 

x = -x-y )-^-x + y/2^n(t) (11) 



p = 2 7 £> 



Jz 



p. (10) 



y = -iv (12) 

(13) 

By solving d^x = for Q, where dn = d/dQ, the feedback which locally mazimizes the 
rate at which the x component grows can be found. This gives the optimal feedback 
strength as a function of x: f2 opt (x) = \J T 2qjx. Using this and the initial condition 
x(0) = the solution to Eq. §\A§ is 



Zopt(t) = Vl - e~^. (14) 
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The optimal feedback strength as a function of time is computed by substituting Eq. ([Til) 



into fi opt : 



This is precisely the evolution of Jacobs' rapid purification protocol gives [12] . Moreover, 
Wiseman and Bouten have shown that this is the globally optimal control protocol for 
this problem [TJ] (see also Ref. [20] ). To calculate the improvement feedback offers 
we solve x m fb(£mfb) = ^opt(^opt) for the ratio of t op tAmfb- Then we take the limit that 
t m fb — > oo to give toptAmfb = 1/5 = 1/2. Thus the asymptotic speed-up (improvement) 
of the ideal feedback over the minimal control (open-loop) protocol to produce a state 
with a given fixed value of x(tf) very close to one is S = 2. Of course, once the desired 
Bloch vector length is reached we may rotate final state to any desired state, as we have 
assumed an effectively instantaneous rotation. Thus we may rapidly prepare any qubit 
state. 



4. Rapid state preparation with imperfections 

We now calculate the effect of four different imperfections or limitations on the rapid 
preparation of an x eignenstate by feedback. 



4-1- Constant feedback strength 

First we consider the case where the feedback strength is time independent (constant). 
It is unclear whether such a strategy can rapidly purify at all, and, if it can, whether 
it affords the same asymptotic advantage as the optimal strategy. In addition this 
analysis is important experimentally for two reasons. Firstly, a fixed feedback strength 
is simpler to implement experimentally. Secondly, the optimal feedback strength (15) 
is unbounded at time zero, so with the analysis here we can probe the usefulness of 
bounded strength controls. 

Based on the asymptotic value of the optimal feedback strength Q opt (oo) = y/2^y, 
we parameterize the constant feedback strength as Q = y/2Fja, where a E (—00,00). 
The x Bloch component becomes 

x — — ary — x + y 27SZ 

= — ary — ja 2 x + 27a. (16) 
Solving this equation with the initial condition x(0) = gives 

x(t) = -^(l-e-^ 2 )). (17) 
1 + cr 

In Fig. [JJ Eq. (j!7p is plotted for different values of a. Negative values of a result in 
the —x eigenstate being prepared. When \a\ > 1 the x component of the Bloch vector 
rapidly increases for t 7 -1 but reaches a steady state value less than one. When 
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\a\ < 1 the Bloch vector increases more slowly, and also asymptotes to a value less than 
one. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



Time (y 1 ) 

Figure 1. The x component of the Bloch vector for different constant feedback 
strategies. The dotted line (blue) is the ideal feedback i.e. Eq. (fT4|) . The thick solid 
line (black) is for the open-loop control protocol, i.e. Eq. ([7]). The remaining lines are 
for the constant feedback strength strategies. The feedback strength is parameterized 
by f2 = ^f2Fja where a £ (—00,00). From the uppermost curve the solid lines (green) 
are for a = 1,3,5 and 15. The dashed lines (green), starting from the upper most 
curve, are for a = 0.5,0.25,0.1 and 0.025. 

The steady state at long-times is x ss = 2a(l + a 2 ) -1 , which is maximized for a = 1, 
in which case 

x (t) = 1 - e" 27 *. (18) 

In Fig. [I] we analyse the short to medium time regime. The optimal constant feedback 
begins to out-perform the no feedback case for t > 0.7687~ 1 . Figure [2] shows the 
long time regime, in which it is appears that Eq. ( TT8|) approaches the the optimal time 
dependent feedback curve ( fT4l) . Interestingly, when a = 0.9 the feedback performs 
worse than the open-loop control until t > 1.537~ 1 . For 1.537" 1 < t < 3.657 -1 feedback 
performs better, and after this interval the feedback again performs worse than the open- 
loop control. This suggests that an optimal state preparation (or purification) protocol 
for bounded strength control (or bang-bang control) might switch between periods of 
feedback and measurement multiple times. 

The only question that remains is whether the optimal constant feedback strength 
protocol has the same asymptotic advantage as the optimal time dependent feedback. 
This can be determined by comparing the ratio of the time taken to a fixed Bloch vector 
length 1— e. Taylor-expanding Eq. (TH|) at long-times gives x opt (t) = 1 — |e~ 27topt . Setting 
this equal to 1 — e and solving for t opt gives the time taken to a fixed Bloch vector length. 
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Figure 2. A comparison of the time it takes to prepare an x eigenstate for four different 
strategies. The dotted line (blue) is for the optimal time dependent feedback - that 
is, Eq. (Tl4| . The thick line (black) is the Bloch vector length for the open- loop control 
protocol. The dot-dash (green) line is the asympototicly optimal constant feedback 
(q = 1), Eq. (fT8|). The dot-dot-dash (red) line is sub-optimal constant feedback - that 
is, Eq. (|17p with a = 0.9. Inset: an exploded version of the same plot. 



This process must be repeated for Eq. ( [18]) ; the time at which Eq 
denoted by t c . The ratio of the times is 
t c lne 
£ opt In e + In 2 
which approaches 1 as e approaches zero. 

4-2. Calibration errors 

Consider now the case in which the applied control has calibration errors. For the 
optimal time-dependent control a time-dependent calibration error could be modelled by 
Q{t) = {■sffi/s/l - e- 2 ^)(l + 5), where 5 G [—1, 1]. Because these errors are systematic 
it is sensible to assume that 5 is constant. Substituting Q(t) into the Bloch equations 
gives 

x = -jx + y/2rfl(t)(\ + 5)- Ul{tf{\ + S) 2 x. (20) 

This has the solution 

x{t) = -e- 2 *T)(l -5 2 ). (21) 

At long-times this asymptotes to x ss ~ 1 — 5 2 /2 for small 5. In Fig. [3] we plot Eq. ( 121]) for 
different calibration error values. A 25% calibration error performs worse than the open- 
loop protocol for t > 2.157 -1 . However, provided the calibration error is less than 5% 
feedback control will out-perform open- loop control for t < 5. For comparison we also 
plot the performance of the asymptotically optimal constant feedback protocol, which 
outperforms the time-dependent feedback with a 5% calibration error for t > 3.057 -1 . 



181) reaches 1 — e is 



(19) 
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Figure 3. The effect of timc-dependent calibration errors on RSP. As before the 
dotted line (blue) is the optimal control; the dot-dashed line (green) is the optimal 
constant (a = 1) control; and the thick line (black) is the open- loop control strategy. 
The three dashed (magenta) lines from top to bottom are for 5 = [0.25, 0.5, 0.75]. The 
dashed line on the inset figure depicts 5 = 0.05. 



4-3. Efficiency and noise on control errors 

Now consider the case in which the experiment is constrained by a detection inefficiency 
n. We analyse this scenario in two ways. First we examine the situation when the 
experimenter applies the control unaware of detection inefficiencies. The second scenario 
is more realistic. Here we assume that there is a detection inefficiency and ask what the 
optimal control in this situation is. Note that detection inefficiency also describes the 
case when there is random noise on the control pi]. We note that the effect of detection 
inefficiency on a related rapid purification protocol was studied in Ref. |21j . 



4-3.1. Oblivious inefficient detection When the system is constrained by inefficient 
detection the equation of motion for the x Bloch component is 

x = -X7 + V2^n - (22) 

If one were oblivious to the detection inefficiency, one would apply the optimal control, 
f2 opt (t), from Sec. EJ The solution of the equation of motion under such a control strategy 
is 



x (t) = ^?R—Ly/i - e-27*. (23) 

Clearly the protocol only works for r) G (0.5, 1] and asymptotically the greatest 
achievable Bloch vector length is x ss = y/2r) — 



4-3.2. Optimal control for inefficient detection To optimize the control for inefficient 
detection one simply takes the derivative of Eq. (122]) with respect to Q and solves to 
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find the optimal control 

«,(«) = (24) 
x 

The solution of the x component using this strategy is 

x(t) = y/tjy/il-e-**). (25) 

Asymptotically the steady state value of this stratergy is a; ss = y/rj. In Fig. H]the steady 
state values of Eq. ( 123]) and Eq. (|25|) . are plotted for all values of r\. The control optimized 
for detection inefficiencies, Eq. (125]) . obviously outperforms the oblivious control case 
for all rj. This is evident in Fig. [51 where Eq. ( 1231) and Eq. ( 1251) are plotted for ^ = 0.85, 
although the two strategies are barely distinguishable. With this efficiency, the state 
preparation becomes worse than that of constant feedback strength, and that of the 
open- loop strategy, for t > 7 -1 . 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

T| (dimensionless) 



Figure 4. A comparison of the steady state inefficiencies for the oblivious control 
scheme, plotted as the solid line (red), and the optimal control for inefficient detection 
dotted line (blue). 



4-4- Time delay 

While Eq. fl2J) is a very general FBME, it assumes instantaneous feedback (that is, 
zero delay time in the feedback loop). It is possible to include the effects of a small 
finite time delay in the analysis while still using a master equation approach [22] . Note 
however that the resulting equation is not guaranteed to be a valid (Lindblad-form) 
master equation and will almost certainly give nonphysical results for r too large with 
respect to 7 _1 . The equation describing the feedback process including a small time 
delay r is [22] 

p = C\p\- iryf^i [F, C [Xp + P X] - X(C [p]) - (C [p])X] . (26) 

Recall that C [.] was defined in Eq. (j2J). Substituting X = J z and F = Q(t)J y into 
Eq. fT26|) . and setting r\ — 1 gives 

p = C\p\- irU{t)y/¥f[J y , £ [J zP + P J Z ] - J Z (C \p\) - (£ \p\)J z \ . (27) 
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Time (y 1 ) 

Figure 5. A comparison of the dynamics of the x component of the Bloch vector 
at relatively short times for the oblivious control scheme (the lower dot-dot-dashed 
line) and the optimal control for inefficient detection (the dashed line); r\ — 0.85. The 
dotted line is the optimal control; the dot-dashed line is the optimal constant (a = i) 
control; and the thick line is the open-loop control strategy. 



From this equation it is straight-forward to calculate the dynamical equations for the 
Bloch components: 

±= - x 7 + y/^yQ - - ^— J: ( 28 ) 



y= - iy (29) 

i = - ^ - 2zVL 2 1 t. (30) 

The feedback strength, Q, must now be chosen in order to rapidly prepare an x 
eigenstate. 

As for the case of inefficiencies, we choose two scenarios to probe time delays in 
RSP. The first scenario, presented in Sec. I4.4.1I examines the performance of a control 
protocol which has time delays, but the experimentalist is oblivious to these delays. The 
second situation the experimentalist is aware of the feedback delay and compensates for 
it, which is presented in Sec. I4.4.2I In both cases we take the feedback strength Q 
to be constant, equal to its asymptotically optimal value. This is necessary to obtain 
analytical results, as obtained in all of the other scenarios in this paper. 

4-4-1 . Oblivious time delay Consider an experimentalist who is unaware of time delays 
in their feedback loop. Due to the cost and difficulties of applying time dependent 
controls the experimentalist decides to apply asymptotically optimal time independent 
control Q = in their RSP procedure. The solution of Eq. ( 1281) in this scenario is 

x(t) = (1 - e" 2 ^) (1 - yr) - (31) 

The steady state solution is x ss = 1 — jr. In Fig. |6]Eq. ( )3T|) is plotted for different time 
delays. 
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Figure 6. A comparison of different feedback strategics for different time delays. 
The dotted line is the optimal control; the dot-dash line is the optimal constant 
(a = 1) control; and the think line is the measurement alone strategy. The dashed 
lines (magenta) from top to bottom are for feedback with oblivious time delay with 
t = 5 x 10 _2 7 _1 , 10 _1 7 _1 . The dot-dot-dashed lines (cyan) from top to bottom are 
for the asymptotically optimal control with time delays r = 5 x 10 _2 7 _1 , 10 _1 7 _1 . 



4-4-2. Asymptotically optimal control for time delays Consider first a local-in-time 
optimization of Eq. (12 8j) . This is again achieved by solving dnx = for Q, which gives 
^opt = — x =t a/x 2 + 127?7r/(3v / 27r). A second derivative test shows that Q^ pt is the 
local maximum, thus we define f2 opt = fl^pt- As Eq. ( !26l) is a perturtabative correction 
to order r, it is only sensible to consider the terms to order r: 

X X 6 



If we restrict, as in the oblivious case, to constant feedbackj§| then optimal 
asymptotic value in the presence of time delays is thus Q = ^/2rf(l — S^t). Under 
this control statergy the solution of x is 

x f t ) = 2{l ~ 37T ) ( 1 _ e -2jt + ^rt-9^t\ (33) 

W 2(1 - 37r) + 97 2 r 2 V J K J 

It is apparent from Fig. |6] that the asymptotically optimal control significantly 
outperforms the oblivious senerio for t > 7~ x . In fact it gives a final Bloch vector 
length of almost unity. This is because the coefficient Eq. ( |33l) . to first order in r, is 
unity. However the terms in the exponent slow the purification, so that no improvement 
over the open-loop case is found at any time even for r = 0.057" 1 . 



5. Comparison of imperfections 

In order to compare the imperfections we determine the allowed parameter ranges for 
the different imperfections that will enable the feedback control case to perform better 

§ The alternative would be to solve the problem numerically. The relevant equation with locally 
optimal feedback (|32[) is x = —2:7 + 777/.T — 2^ 2 ij 1 t / 'x 3 . 
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than the minimal control (open-loop) protocol at t — 2^^ 1 and t = 47 -1 . Although 
the choice of times is arbitrary, the first time chosen is, roughly, when the long-time 
limit (t ^> 7 _1 ) approximation begins to be reasonable, while the second is well in the 
asymptotic regime. Table [TJ summarizes the results. 



Control Imperfection 


Parameter Range A 


Parameter Range B 


Constant FB strength 


0.956 < a < 1.189 


0.9415 < a < 1.0678 


Time dep. cal. errors 


< 5 < 0.165 


< 5 < 0.0659 


Inefficient det. oblivious FB 


0.973 < 77 < 1 


0.9956 < 77 < 1 


Inefficient det. optimal FB 


0.972 < 77 < 1 


0.9956 < 77 < 1 


Time delay (7 _1 ) oblivious FB 


< r < 0.0045 


< r < 0.0020 


Time delay (7" 1 ) optimal FB 


< t < 0.0146 


< r < 0.0195 



Table 1. The parameter ranges for which imperfect feedback performs better than the 
minimal control protocol at a fixed time. The first column describes the type of control 
imperfection. The second column, Parameter Range A, gives the approximate values 
for the parameter (imperfection) under consideration at t = 27 _1 . The third column, 
Parameter Range B, gives the approximate values for the parameter (imperfection) 
under consideration at t = 47 



Control Imperfection 


Parameter Range A 


Parameter Range B 


Constant FB strength 


0.922 < a < 1.135 


0.9860 < a < 1.0141 


Time dep. cal. errors 


< 5 < 0.1246 


< 5 < 0.01412 


Inefficient det. oblivious FB 


0.9846 < 77 < 1 


0.9998 < r/ < 1 


Inefficient det. optimal FB 


0.9844 < 77 < 1 


0.9998 < 77 < 1 


Time delay (7" 1 ) oblivious FB 


< r < 0.00547 


< t < 0.0000996 


Time delay (7 -1 ) optimal FB 


< r < 0.02598 


< r < 0.004612 



Table 2. The parameter ranges for which imperfect feedback performs better than 
the minimal control protocol at a fixed Bloch vector length. The first column 
describes the type of control imperfection. The second column, Parameter Range 
A, gives the approximate values for the parameter (imperfection) under consideration 
at x = 1 — 1CP 2 . The third column, Parameter Range B, gives the approximate values 
for the parameter (imperfection) under consideration at x = 1 — 1CP 4 . 

From Fig. [TJ it was clear that a ^> 1 is required for short time evolution to perform 
like the open-loop evolution, but the dynamics induced by a large a caused the steady 
state Bloch vector length to be poor. Thus the second row in Table [TJ suggests that 
strong feedback is not necessary as a is of order 1. This implies that bounded strength 
control protocols can work in RSP. The third row in Table (OQ) suggests that the feedback 
protocol is not very sensitive to calibration errors, as a 16.5% error in the calibration 
seems rather high. 
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The fourth and fifth rows of Table [T] show the effect of detection inefficiencies 
on RSP. Even in quantum optics with high-efficiency homodyne detectors, the overall 
efficiency is seldom greater than r\ < 0.9 (see e.g. Ref. [23]). Thus, unfortunately, the 
RSP protocol is likely to be severely compromised by detection inefficiencies, regardless 
of whether one is aware of them or not. These rows also represent the effects of white 
noise on the control. 

The final row indicates that RSP is most sensitive to time delays. Outside the 
small t regime it is unlikely that this sort of RSP will be feasible using Markovian or 
Bayesian feedback]]]] If, however, a particular physical system has 7 _1 which is of the 
order of a Meghertz (or smaller) then feedback delays due to classical circuitry could be 
made much smaller than 7 _1 . In this case RSP is remarkably robust to time delays, in 
the sense that the asymptotic x ss differs from unity only in second order in the delay r. 

In Table [2] we consider the parameter ranges that allow feedback to perform better 
than the minimal control protocol at a fixed Bloch vector length. We consider two 
lengths: x — 1 — 10~ 2 (the second column) and x = 1 — 10~ 4 (the third column). The 
same trends described above also hold for this table. 

6. Discussion 

The most important conclusion of our analysis is that quantum feedback control for 
rapid state preparation is quite sensitive to detection inefficiencies. In particular, locally 
optimizing the feedback, taking into account the efficiency 7], limits the length of the 
Bloch vector in the long-time limit to ^Jrj. However, our analysis also turned up some 
encouraging conclusions. Surprisingly, control in an unbiased basis seems robust to 
delays in the feedback loop, at least in the sense that the asymptotic Bloch differs from 
unity only in second order in the delay r. 

The methodology presented in this article, approximating Bayesian feedback with 
Markovian feedback, is not restricted to studying imperfections in the control. It also 
allows one to study the performance of control protocols with system imperfections. For 
example, consider performing RSP when the system has isotropic dephasing noise at 
rate r iso and the decay rate from the excited state is T d . The FBME in this situation 
is: 

p = c [ P ] + T ™ v N p + T * v p- ( 34 ) 

From this equation it is easy to derive an equation for x(t) 

x (t) = J 2 ^ V — y/i _ e (-2t 7 -8fr iso -r d ^ (35) 

w V 27 + sr iso + r d v 1 ; 

Here we are modelling locally optimal rapid state preparation, with detection inefficiency 
included as a control imperfection. The main point is that, by approximating Bayesian 

| We note, however, that the open-loop scheme presented in Ref. 0] allows sub-optimal RSP with a 
single conditional unitary at the end, as in the minimal control protocol of Sec. [21 



Quantum feedback for rapid state preparation in the presence of control imperfectionslh 

feedback with Markovian feedback we can begin to analyse analytically the performance 
of feedback protocols with many different control or system imperfections. We note that 
the effect of x dephasing and decay from the excited state was studied in Ref. |24|.The 
advantage of our method is two fold: firstly our analysis is algebraically simpler; and 
secondly we derive time-dependent solutions, not just steady state values. 

In future work we think it would also be interesting to see if the Markovian approach 
to rapid x-eigenstate preparation presented here could be generalized to all qbit states 
(by analogy with the stabilization of a qbit under monitored spontaneous emission in 
Refs. [251 IE]) with and without imperfections. 
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